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The Watts-Strogatz algorithm transferring a regular lattice to the small world network is modified 
by introducing preferential rewiring constrained by connectivity demand. The probability to link 
to/ unlink form a node is dependent on a vertex degree and adjusted by some threshold. For 
each threshold value there exists a probability at which the resulting stationary network has degree 
distribution with power-law decay in large interval of degrees. 

PACS numbers: 87.23.Ge, 89.75.Hc, 89.20 Hh 



I. INTRODUCTION 

Networks are used extensively to study and describe 
topology in various complex real systems. From commu- 
nication networks like World Wide Web, through nets 
of social relations like network of friends or collaborat- 
ing scientists, to biological systems like protein network, 
neural network or cell metabolism — a similar structural 
organization has been discovered, see Q, 0, Q for review 
of data analysis and bibliography. 

The simplest and the most intensely studied charac- 
teristic of a network is degree of a vertex. Degree, k, of 
a vertex is the total number of its connections. A physi- 
cist would say degree is actually the number of nearest 
neighbors of a vertex. Distribution of vertex degrees of an 
entire network is the network's basic statistical character- 
istic. Although the degree of a vertex is a local quantity, 
the degree distribution often determines some important 
global features. For example, in case of a regular lattice 
(a network explored in many physical applications) the 
degree distribution is concentrated at a single value what 
exactly proves the homogeneity of the network topology. 

The studies of other topologies began with the random 
graph theory of Erdos and Renyi |4| . The proposition of 
Watts and Strogatz |5j that followed, called small world 
network, captures best the features of both regular lattice 
and random graph. The algorithm may be summarized 
as follows: begin with a regular lattice, then rewire an 
edge that is chosen with some probability p. The crucial 
properties of the resulting network are that it is highly 
clustered like a regular lattice, and at the same time it 
has small path lengths like a random graph. 

The topology of some graphs found in technology such 
as electronic circuits, where electronic components are 
represented as vertices and the wires are edges, can be 
viewed as a class of such networks The magnetic 
materials, including Ising models, based on small world 
networks are investigated widely in the aim of construct- 
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ing magnetic nanomaterials |7j, |8( . The neural networks 
may also be considered as kind electronic circuits 0, 
according to the idea of 'save wiring' as an organizing 
principle of the brain [Tl| . 

The degree distribution of small world network is sim- 
ilar to a random graph. It has pronounced peak at the 
average degree < k > and decays exponentially for large 
k. Thus the topology of the network is relatively homoge- 
neous. However, the degree distribution of most real net- 
works significantly deviates from such distribution. For 
large number of real networks degree distribution has a 
power-law form: 

P{k) oc kr' 1 (l) 

where 7 is the degree exponent. A network with the 
property Q is called scale free network. In a scale free 
network, the smaller value of 7, the greater role of, so- 
called hubs: vertices with large degree. Especially, if 7 < 
3 then the surprising network features such as robustness 
against accidental failure of a vertex occur ^1 • 

Barabasi and Albert proposed the method for con- 
structing a network with the power-law degree distri- 
bution [l3j. The two basic features of the model arc: 
the network growth (in each time step a new vertex is 
attached to a graph), and the preferences in the attach- 
ment (a vertex is in favor according to its degree). Doro- 
govtsev and Mandes 01 proved that aging (a preference 
growth based on the age of the node) could lead to scale 
free network. Dangalchev implied that the scale free 
topology can be restored in a model with a static number 
of nodes. However, the organizing principles that govern 
and are common in a wide range of complex networks are 
far from being fully understood |16| . 

In the following paper we present the effect to the 
degree distribution of adding preferences to rewiring in 
Watts-Strogatz algorithm. In analogy to the small world 
algorithm, we start from a regular lattice (a square lat- 
tice) and set the preferences in rewiring to be related to 
a vertex degree and adjusted by some threshold. In addi- 
tion, we set the constraint of preserving the connectivity 
of a graph. The resulting network at properly adjusted 
parameters has the degree distribution of power-law type 
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in the wide interval. The degree exponent 7 is less or 
close to 2. Then in a tail, the distribution becomes an 
exponential one. 

The algorithm is presented in Section 2. Section 3 con- 
tains results from simulations. The planned development 
of the model is proposed in Section 4. 

II. THE ALGORITHM 

The vertices of the lattice are numbered by: 
1,2,..., iV. Each edge is characterized by the two num- 
bers from and to representing two vertices linked by this 
edge. The graph is represented as the vector of size N of 
lists of vertices — neighbors of subsequent vertices. The 
graph is not considered as directed, though it can easily 
be adapted to a directed one. 

There are two basic parameters in the model: 

— probability to rewire an edge: p 

- threshold to set in favor a vertex in the preference 
function: T 

At each time step, N vertices are considered for 
rewiring. A vertex for rewiring is selected at random. 
The basic procedure, called, EdgeEvolution, requires 
three parameters: from, p and T. The procedure in- 
volves considering changes in all edges linked to a vertex 
from. The rewiring an edge means change a present 
value to to some other. The steps in the procedure are 
described below; the illustration of the 'before and after' 
is shown in Fig. ^ 

EdgeEvolution ( from, p, T) : 

(a) for each vertex to from the list of edges of the vertex 
from do 

(b) if degree to is greater than 1 

(c) then 

(d) let £1 £ [0, 1] be a random number 

(e) if 



degree of to 

(f) then accept vertex to for unlinking 

(g) for each vertex to accepted for unlinking 

(h) let newdo £ {1, 2, . . . , N} be a random number 
but new-to ^ from 

(i) let £ 2 € [0, 1] be a random number 
0) if 

degree of newdo 

£2 < — y 

(k) then accept newdo: 

edge(from,to) := edge(from,new_to) 




before after 



FIG. 1: An edge linking vertices from and to ( before graph) 
is rewired to the edge between vertices from and newJto (af- 
ter graph) 



(1) otherwise go to (h) 



Remarks: 

— The algorithm conserves both the number of vertices 
and the number of edges. Steps (a)-(f) prepare the list of 
edges of the vertex from to be rewired, while steps (g)- 
(1) fix new connections. Each edge accepted for rewiring 
has to be rewired (see the loop(l)). 

— The preferences (e) and (j) are chosen to prevent un- 
linking and favor linking to vertices with degrees greater 
than the threshold value. 

— The condition (b) is crucial to conserve connectivity 
of a graph. Any vertex with the degree equal to 1 cannot 
be unlinked. Otherwise we face the problem of dramati- 
cally increasing number of isolated nodes. 

— A randomly selected new vertex for linking to must 
be different from the number from to avoid loops in a 
graph, line (h). 

The update of information regarding the vertex degree 
is made after each time step. Therefore, if the network 
is evolving with small p, then the evolution may be seen 
as asynchronous, while with large p, the evolution may 
look as if performing synchronously. 



III. RESULTS 



We test by simulations the algorithm on a square lat- 
tice with L = 100 and L = 200. It means that number 
of vertices considered is N = 10 4 ,4 * 10 4 , respectively. 
The number of edges is AN. We observe no difference in 
results obtained for different lattice sizes. 

For any p and T in a few time steps the initial 5- 
like degree distribution centered at k — 4 spreads to 
some peaked distribution with the center around 4. As 
our focus is on stationary distributions we give the sys- 
tems time to achieve stable states. If preferences are not 
switched on then the network quickly, in less than 100 
steps, reaches the stationary state with the exponential 
degree distribution . If the preferences are switched on 
then the system requires much more time to stabilize. 

The degree distribution is a two-part function in the 
stationary state. The first part represents properties of 
vertices which are not preferred, i.e., vertices with de- 
gree smaller than the threshold (k < T). The shape of 
this part can be approximated by a power-law with the 
degree exponent 7 greater than 3. The second part of 
the distribution represents tail properties, i.e., probabili- 
ties to meet vertices with degrees much greater than the 
threshold T, k >> T. The tail decays exponentially. Be- 
tween these two separated regions there is a transitional 
region which can also be characterized by the power-law 
dependence on vertex degree. It appears that for the 
given T there exists a certain value of p (denoted pr) 
for which the first region and the transitional region be- 
have according to the common power-law decay with the 
exponent much smaller than 3. The distribution at pt 
is critical — for p > pt, the parts of the distribution 
become indistinguishable. 

On Fig. |21 we present the degree distributions for var- 
ious values of p and T = 8. We observe the evolutions 
up to 50000 steps and find that after 10000 steps the 
systems are in stationary states — the degree distribu- 
tions are identical. The distributions are shown on log- 
log plots for a better visualization of their power-law 
features. The regions discussed above are easily iden- 
tified. For the degrees smaller than the threshold, here 
k < 8, the probability decays in a polynomial way. For 
the degrees large — tail degrees (fc > 50), the probability 
decay is exponential. The transition region is between 
these two parts, i.e., for 8 < k < 50. When p rs 0.90 
the widest region with the power-law dependence is ob- 
served, 2 < k < 50. Therefore, p% w 0.90. Notice, that at 
p = 0.1 there are vertices with the highest degrees of all 
presented results in Fig. |2 The network is highly inho- 
mogeneous. However, the highest degrees observed by us 
at this threshold value occur at p = 0.30 : k max « 500. 
The transition region with this p characterizes by degree 
exponent 7 = 0.43(0.97). 

For comparison on Fig. [3] we show the critical degree 
distributions for T = 12 and T = 4 respectively. We 
found pi ~ 0.7 andpi2 ~ 0.9. Generally, we observe that 
the greater value of T, the higher degrees are attainable. 
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FIG. 2: On subsequent figures the log-log plots of the degree 
distributions with T = 8 and different p are presented. The 
degree exponent values are calculated for points marked by 
gray big dots and are written at the plots together with their 
Pearson coefficients r 2 . 



In case of T = 12 and p — 0.05 vertices with a degree 
greater than 600 happens. If T = 16 then vertices with 
k = 1000 are present on a network, see Fig. Notice, 
that in this case the transition region spreads from k = 16 
to almost k — 200. It implies that with further growth 
of the threshold the two-part distribution splits into two 
separate curves — there is no transition region. 



IV. MODEL DEVELOPMENT 

Watts and Strogatz showed that small worlds captures 
the best of two graph worlds: the regular and the ran- 
dom. However, there has been no work discussing effect 
of preferences in rewiring. Our work is a step into this 
direction. In particular we ask whether small worlds can 
lead to scale free networks. In the limit of long runs 
(many time steps) the evolving network with preferential 
rewiring may be self-organized into free scale structures. 
The driving feature is the requirement to preserve net- 
work connected. 

Our work supports the idea that not growing network 
may underlie the formation of scale free network. In fu- 
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FIG. 3: Log-log plots of critical degree distributions for T — 
12 and T = 4. 




FIG. 4: Log-log plot of the widest degree distribution ob- 
tained in simulations with T = 16. 



ture attempts to gain more understanding of the orga- 
nizing principle in complex systems we will study other 
preference functions. In particular, we will investigate 
the modification of the rewiring procedure that restricts 
the rewiring to a local rule — the set of accessible ver- 
tices when rewiring a (from, io)-edge will be restricted to 
the set consisting of nearest-neighbors of a to- vertex. By 
this restriction we will be able to consider new features 
arising in cellular automata systems due to the network 
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evolution. Cellular automata are the simplest systems 1472/PO3/2003/25 
used widely in modeling different natural systems, see 

m 
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